Cooperative optical non-linearity due to dipolar interactions in an ultra-cold Rydberg ensemble 
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We demonstrate a cooperative optical non-linearity caused by dipolar interactions between Rydberg atoms in 
an ultra-cold atomic ensemble. By coupling a probe transition to the Rydberg state we map the strong dipole- 
dipole interactions between Rydberg pairs onto the optical field. We characterize the non-linearity as a function 
of electric field and density, and demonstrate the enhancement of the optical non-linearity due to cooperativity. 
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Photons are robust carriers of quantum information and 
consequently there is considerable interest in the development 
of photonic quantum technologies. As optical non-linearities 
are extremely small at the single photon level H) attention 
has focussed on linear optical quantum computing ||2][3). In 
parallel, work has been carried out on materials with a large 
Kerr effect JU [5] |6] [7] |8) potentially enabling non-linear pho- 
tonic devices. Theoretical work has explored some of the dif- 
ficulties in realizing a high fidelity quantum gate based on the 
Kerr effect [9j. An alternative mechanism for generating an 
optical non-linearity, for example a cooperative non-linearity 
due to dipolar interactions, could open new avenues for pho- 
tonic quantum gates [ 1 1 . In a dipolar system the electric field 
is modified due to the local field of the neighbouring dipoles 
ifTTl . Such local field effects can give rise to cooperative be- 
haviour such as superradiance lfT2l \l3l and optical bistability 
HUE]. 

In this paper we demonstrate a cooperative optical non- 
linearity due to dipole-dipole interactions between Rydberg 
atoms. These strong interatomic interactions are sufficient 
to prevent excitation of neighbouring atoms to the Rydberg 
state lfl6l . This gives rise to a blockade mechanism which 
has been observed for a pair of trapped atoms IfTTl [T8l and 
an atomic ensemble [19|. In our work the effect of strong 
interactions between Rydberg pairs is mapped onto an opti- 
cal transition using electromagnetically induced transparency 
(EIT) [20 21 J. The resonant dark state responsible for EIT is 
modified by the dipole-dipole interactions, causing suppres- 
sion of the transparency on resonance. The resulting optical 
non-linearity depends on interactions between pairs of atoms 
and is a cooperative effect where the optical response of a sin- 
gle atom is modified by the presence of its neighbours. 

To show how dipole-dipole interactions give rise to a co- 
operative non-linear effect, we consider the atom pair model 
ll22l shown in fig. 1(a) for three level atoms with ground \g), 
excited |e), and Rydberg |r) states. These states are coupled 
by a probe laser with Rabi frequency f2 p and a strong coupling 
laser with Rabi frequency il c . In the non-interacting case with 
probe and coupling lasers tuned to resonance the dark state is 
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FIG. 1. (Color online), (a) Pair model for interacting three-level 
atoms. States \g), \e) and \r) are coupled by a probe field Q p and a 
strong coupling field Sl c . Dipole-dipole interactions couple the Ry- 
dberg states with strength V(R). (b) Schematic of experiment. EIT 
spectroscopy is performed on an ultracold 87 Rb atom cloud. 



where tan 8 — £! p /f2 c and cf> r is the relative phase between 
probe and coupling lasers. This state is not coupled to the 
probe field, leading to 100 % transparency independent of the 
mixing angle, 8. 

Dipole-dipole interactions modify this picture. The effect 
of interactions is to detune the doubly excited Rydberg state 
\rr) by the energy V(R), which depends on the interatomic 
separation, R. In the limit of strong interactions (V > il c ) the 
detuning becomes sufficient to block excitation of \rr). This 
means \iPd) is no longer a valid eigenstate, with the new zero 
energy eigenstate given by [24]: 
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This is no longer a dark state as the interaction has mixed in 
the rapidly decaying |ee) state in place of \rr). The result is 
that the EIT is suppressed by an amount which depends on 
the probe electric field through 8. The change in the transmis- 
sion, T, of the probe laser is due to a change in the imaginary 
part of the complex susceptibility at the probe laser frequency, 
Xi = Im[x(w p )]. These are related by T = 1 — exp{— k\i£\ 
where k is the wave-vector and t the path length through 
the medium. The transition between the non-interacting state 
eq. ([T| and the blockaded state eq. |2} depends on the dipole- 
dipole coupling strength V. As a result, the single atom sus- 
ceptibility depends on both probe electric field, giving an op- 
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FIG. 2. (Color online). Suppression of EIT due to dipolar interac- 
tions for |r) = 48S 1/2 at TV = 2.2 x 10 10 cm" 3 . Thin lines show 
absorption data at fi p /27r = 0.3 MHz (a) and 1.0 MHz (b), revealing a 
suppression in transparency on resonance. Theory curves calculated 
from the pair model are plotted for V/2n — 10 MHz (thick line) and 
V/2-k = (dashed line) using 48Si/ 2 state lifetime of 58 ^s (28), 
tt c /2-K = 2 MHz, 300 kHz probe laser linewidth and 150 kHz relative 
linewidth between probe and coupling lasers. 

tical non-linearity, and also the distance to neighboring atoms, 
which gives a density dependent cooperative non-linearity. 
Hence the absorption and refractive index contributed by any 
atom depends on its proximity to other atoms, leading to a 
non-linear density scaling. 

Our experiments are performed on a laser cooled 87 Rb 
atom cloud using the experimental setup described in |25| 
and shown schematically in fig. 1(b). Atoms are loaded into 
a magneto-optical trap (MOT) for 1 s after which the cool- 
ing light is extinguished. Atoms are then prepared in the 5s 
2 Si/2 \F — 2, nip = 2) state (\g)) by optical pumping on 
the 5s 2 S 1/2 F = 2 — > 5p 2 P 3/2 F' = 2 transition. By vary- 
ing the amount of optical pumping, the fraction of atoms in 
F = 2 and hence the density in \g) can be controlled with- 
out changing the cloud size. Counter-propagating probe and 
coupling lasers are turned on simultaneously to perform EIT 
spectroscopy as shown in fig. 1(b). The probe beam absorp- 
tion is detected using a photodiode with a 20 MHz bandwidth. 
The coupling laser has a power of 90 ±5 mW and is frequency 
stabilized on resonance to the 5p 2 P 3 / 2 F' = 3 — » ns 2 Si/2 
transition using an EIT locking scheme [26|. The probe laser 
drives the 5s 2 Si/ 2 F = 2 — > 5p 2 P 3 / 2 F' = 3 transition 
and is scanned across the resonance from A p /27r = —20 — > 
+20 — * —20 MHz in 960 /is. This double scan technique 
provides useful information on atom loss |25|. The probe and 
coupling lasers are circularly polarized to drive a + -a~ tran- 
sitions to maximize the transition amplitude to the Rydberg 
state. The path length through the ensemble along the probe 
axis was measured by fluorescence imaging of the cloud after 
preparation in F — 2, giving £ = 0.52 ± 0.03 mm. Using 
the path length, the transmission is converted to the imagi- 
nary part of the susceptibility \i from the relation above. To 
measure the density of atoms in the probe beam, transmission 
data is recorded with the coupling laser off. This is then fit- 



ted using the analytic absorption profile [27], giving a peak 
density of 7V = 2.2 ± 0.2 x 10 10 cm" 3 . At the peak 
density achieved in the experiment the average pair separa- 
tion is R ~ 2 /jm. For states with principal quantum number 
n < 50, this separation corresponds to van der Waals interac- 
tions V(R) = —Cq/R 6 (29). In this regime interactions are 
dominated by nearest neighbours, closely approximating the 
pair model in fig. 1(a), unlike for resonant dipole potentials 
oc 1/i? 3 for which all atoms contribute [ 30 1 . 

Figure 2 shows individual absorption spectra recorded for 
|r) = 48S]y 2 at two probe powers for the first scan through 
resonance. This shows a significant drop in transparency as 
the probe Rabi frequency is increased, demonstrating a strong 
suppression of the EIT resonance as expected. 

To ensure this is due to dipole-dipole interactions and not a 
Stark shift of the Rydberg state due to ion creation, we con- 
sider the effect of ionization on the EIT lineshape. A Monte- 
Carlo simulation is used to model the Stark shift experienced 
by atoms due to the electric field of ions within the cloud. 
A random distribution of atoms at a density /V is generated 
and a certain fraction are randomly chosen as ions. The Stark 
shift for each atom is obtained by calculating the local electric 
field due to all ions in the cloud. The simulated EIT profile 
is then calculated by integrating over the Stark shifted line- 
shape of each atom. Figure 3(a) shows the simulated profile 
for \r) = 48Siy 2 calculated using a polarisability of a^= 
38.2 MHz/(V/cm) 2 PP . This shows that to see a suppres- 
sion of similar magnitude to fig. 2(b) due to ions, there would 
also be a significant shift of the resonance (>1 MHz). Figure 
3(b) shows experimental measurements of the two-photon res- 
onance position as a function of probe intensity, where no shift 
is seen for the experimental sequence. As a further check we 
deliberately create ions by reversing the probe scan direction. 
Ions are created on the blue side of the two-photon resonance 
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FIG. 3. (Color online). Effect of ions (a) Theoretical spectra for 
\r) = 48Si/ 2 for ground state density 2.2 x 10 10 cm" 3 for 0, 2 and 
5 % ion density showing a Stark shift in the two-photon resonance 
before significant suppression of the resonance, (b) Measured two- 
photon resonance for experimental scan (■) shows suppression not 
due to ions. Reversing the probe scan direction to first create ions 
shows a visible Stark shift in the second part of the scan (•), showing 
the sensitivity to ions. Each point is the average of 20 measurements. 
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FIG. 4. (Color online). Optical non-linearity measured for 
|r) = 42Si/ 2 at several densities where N = 2.2 x 10 10 cm -3 , 
each point the average of 20 measurements. High density data is fit 
to a cubic non-linearity up to Q p /2tv = 0.7 MHz, shown by dashed 
line. This gives x' 3 ' ~ 2 x 10 -7 m 2 V -2 . Solid lines show theory 
curves integrating pair-model over distribution of nearest neighbours 
for each density, giving good qualitative agreement. 



due to repulsive pair potentials [32], leading to an observable 
Stark shift in the second part of the scan. Since this shift is 
not seen for the forward probe scan, the suppression observed 
in fig. 2 is due to dipole-dipole interactions. 

We now compare the experimental observations to the- 
oretical predictions of the effect of dipole-dipole interac- 
tions using the pair model of fig. 1(a). To reproduce 
the transmission spectra shown in fig. 2, it is necessary 
to calculate the complex susceptibility at the probe fre- 
quency. This was done using the Liouville equation for 
the density matrix of the two atom system using the basis 
states^), \ge), \eg), \gr),\rg), \ee), |er), \re), \rr)}, where 
states are labelled \j)2- The decoherence terms included 
in the model are the spontaneous decay rates from the Ryd- 
berg and excited states as well as the linewidth of the probe 
and coupling lasers. For the Rydberg states the decay rate is 
a few kHz, however the relative linewidth of the probe and 
coupling laser is approximately 150 kHz and dominates the 
linewidth of the EIT feature. 

Dashed lines on fig. 2(a) and (b) show theoretical absorp- 
tion profiles for the non-interacting case (V = 0), calcu- 
lated using parameters optimized to fit the weak probe data 
at low density where interactions can be neglected. In the ab- 
sence of interactions, the transparency on resonance increases 
with probe Rabi frequency due to population transfer from the 
ground state to the Rydberg state, reducing absorption. This 
is the opposite of what is observed experimentally, as the in- 
teractions reverse the sign of the optical non-linearity. Thick 
lines on fig. 2 show the absorption calculated in the blockaded 
regime (V ^S> Q c ) for V/2ir = 10 MHz, where increasing 
V further gives no additional suppression. The model repro- 
duces the enhanced absorption, however it underestimates the 
magnitude of the suppression. For |r) = 48Si /2> the block- 
ade radius — v/CgA^c ~ 4 /jm 11331 corresponding to 



an average of 5 atoms per blockade sphere. This means that 
we cannot only consider isolated pairs of atoms and the pair 
model does not give a complete description in this case. If 
more than two atoms contribute to the blockade then the re- 
sulting suppression increases due to a higher fraction of pop- 
ulation in |e). 

To explore the cooperative nature of the non-linearity, data 
were taken for \r) = 42Si/2 which has a smaller blockade ra- 
dius of 3 /im J33). This gives an average of 2 atoms in each 
blockade sphere, which should lead to better agreement with 
the pair model. EIT spectra are taken as a function of probe 
power for several densities, in each case taking 20 single mea- 
surements and extracting the average susceptibility on reso- 
nance. Figures 4 and 5 show the measured non-linearity, with 
the susceptibility scaled by density to give the average sus- 
ceptibility of a single atom Xi/N. The optical non-linearity 
is shown in fig. 4, showing that there is a rapid increase in 
susceptibility with probe field until it saturates at Xi/N — 
1 .2 x 10 -20 m 3 once the blockade regime is reached. To quan- 
tify the strength of the non-linearity, a third-order susceptibil- 
ity is fitted to the data up to saturation, shown as a dashed line 
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on the figure. The fit gives x^ 3 ' ~ 2x 10" 
in magnitude than slow light in a BEC [4 1 but with two orders 
of magnitude lower density. This comparison illustrates the 
cooperative enhancement effect arising due to strong dipole- 
dipole interactions. 

For a standard optical non-linearity the single atom re- 
sponse Xi/N varies with electric field but is independent of 
density. In the cooperative case the response also changes 
with density. This is illustrated in fig. 5. For a low probe field 
f2 p /2"7r = 0.3 MHz the population remains in the ground state 
and is insensitive to the dipole-dipole coupling of the Rydberg 
state. Increasing the probe power causes the interaction effect 
to become important and the single atom response becomes a 
function of N . 

To compare the measurements to the theoretical non- 




TV (xl0 lu cm- 3 ) 



FIG. 5. (Color online). Cooperative non-linearity measured for |r) = 
42Si/2, showing strong density dependence at high probe field. The 
interactions increase rapidly with N due to the 1/R 6 scaling, causing 
the susceptibility to plateau as V > Q c . Solid lines show model 
integrated over nearest neighbour distribution. 
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linearity predicted by the pair model it is necessary to inte- 
grate over the range of interactions strengths in the ensem- 
ble. This was done by taking the nearest neighbour distribu- 
tion for a uniform density with a Poissonian atom distribution 
p(R,N) = AttNR 2 exp[-|7r7Vi? 3 ]. For each separation R 
the pair model was used to calculate the resonant suscepti- 
bility at a given probe Rabi frequency. We obtained the to- 
tal susceptibility by summing the susceptibility weighted by 
p(R, N) for each separation. This was repeated for each probe 
Rabi frequency to give the solid curves shown on figures 4 and 
5, calculated using the same parameters as before except as 
n=42, fi c /27r = 2.5 MHz and the Rydberg lifetime is 41 fis 

ma. 

Figures 4 and 5 show the pair model agrees well with the 
data in the limit of low density, where interaction effects are 
weak, and large density or strong probe power, where the sup- 
pression reaches saturation due to blockade effects. As the 
saturation value of the single atom susceptibility Xi/N agrees 
well with observations, this implies that modeling the system 
as an ensemble of isolated pairs is a valid description of the 
interactions, which can be used to predict the behaviour of the 
system. For intermediate probe powers the model underesti- 
mates the onset of suppression, seen in fig. 4, as the predicted 
optical non-linearity is much smaller than the measured value. 
This discrepancy could be explained by a modification to the 
nearest neighbour distribution p(R, N) to consider a Gaussian 
rather than uniform density distribution. 

The data presented here show that the cooperativity can be 
controlled using both probe electric field and density. How- 
ever it is also possible to modify the single atom response 
by changing the strength of the coupling field. As the cou- 
pling strength is increased, not only does the transparency at 
low density increase but stronger interactions are required to 
blockade excitation of Rydberg pair states. This means the 
density at which the cooperativity saturates can be increased. 

In summary we have shown how dipole-dipole interactions 
between a pair of Rydberg atoms leads to a cooperative op- 
tical non-linearity. This occurs due to modification of the 
resonant dark state giving enhanced absorption on resonance 
and a modification of the dispersive properties of the medium. 
The cooperativity has been demonstrated through a change in 
the average single atom response with both density and probe 
electric field. On resonance, the effect of interactions is to 
increase the absorption of the medium, which is a dissipa- 
tive process. This could be avoided by detuning off -resonance 
where the interactions induce a predominately dispersive shift. 
In future work, this cooperative non-linearity will be used to 
map the blockaded atomic excitation onto the photon field, 
enabling the generation of non-classical states of light. 
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